Photoluminiscence of a Quantum Dot Hybridized with a Continuum of Extended 

States 
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We calculate the intensity of photon emission from a trion in a single quantum dot, as a function 
of energy and gate voltage, using the impurity Anderson model and variational wave functions. 
Assuming a flat density of conduction states and constant hybridization energy, the results agree 
with the main features observed in recent experiments: non-monotonic dependence of the energy 
on gate voltage, non-Lorentzian line shapes, and a line width that increases near the regions of 
instability of the single electron final state to occupations zero or two. 

PACS numbers: 73.21. La, 78.67.Hc, 75.20.Hr, 72.15.Qm 



Semiconductor quantum dots (QD's) have attracted 
much attention recently and have been proposed for nu- 
merous applications. It has been established that they 
are promising candidates in quantum information pro- 
cessing, with spin coherence times of several microsec- 
onds ju [ij and fast optical initialization and control 
QD's are also of interest in the field of spintron- 
ics (electronics based on spin) [H [jl §]. Optical means 
of manipulation Q and detecting ]t| the spin were pro- 
posed. QD's might also be used as semiconductor opti- 
cal amplifiers . Most of the above cited works involve 
optical transitions which create or destroy excitons. In 
particular, when a positive gate voltage is applied, con- 
figurations with one more electron can be stabilized, and 
the optical transition takes place between a state with 
the valence band full and only one electron in the con- 
duction band, and a state which consists of two conduc- 
tion electrons and a hole in the valence band (a trion) 
0, i, S, EH EE El El- The gate voltage V g allows 
to control the spin dynamics and the optical emission 

0,E1E3,E1. 

An interesting aspect of the photoluminiscence (PL) 
decay of the excitons and trions, is the manifestation of 
the hybridization of the spin localized in the QD with a 
continuum of extended states, and the related Kondo ef- 
fect 12, [l3|, LLJ] . This hybridization can be described by 
the impurity Anderson model [HI [see Eq. ([7]) below] . In 
the limit of small hybridization and odd number of elec- 
tron in the QD, this model reduces to the Kondo model, 
in which the localized spins have an exchange interaction 
with the spin of the extended states 12j. The Kondo 



physics is also clearly present in transport properties of 
QD's 0,0. 

In recent experiments, the PL that results from the 
decay of the X l ~ trion has been measured as a function of 
V g [13| . The PL is a consequence of an optical transition 
from the trion to a state with one localized electron in 
the QD, hybridized with a continuum. As a consequence 
of this hybridization, the PL is broad (half width at half 
maximum w ~ 1 meV) near the low V g limit of stability of 



the trion (for lower V g , the neutral exciton X° is stable). 
This limit is slightly below V g °, where we call Vg the gate 
voltage for which the configurations with n and n + 1 
electrons in the QD, and no holes in the valence band, 
are degenerate. The striking observed behavior can be 
summarized in three main results: i) as AV g = V g — V g ° 
increases from negative values, the PL first blueshifts and 
then redshifts, ii) the PL line shape is non-Lorentzian, 
with a low energy tail, iii) while w *~ 1 meV for V g ~ V g ° 
or V g ~ Vg 1 (near the limits of stability of one electron in 
the QD), the PL line is much narrower for intermediate 
values of V g . This has been noted before [3|: while at 
intermediate gate voltages Atature et al. succeded to 
prepare a spin state with more than 99.8 % fidelity, the 
coupling of the spin with the reservoir at the edges (V g ~ 
Vg or V g ~ Vg 1 ) spoils this efficiency. 

The line shape has been quantitatively fit only for 
Vg ~ Vg using a one-particle approach with additional 
questionable assumptions [13j |. In this Letter, we calcu- 
late the PL on the whole range of gate voltages using 
many-body variational wave functions (VWF's) 15, lsj . 
Our results provide an explanation of the three main ex- 
perimental results mentioned above and show that they 
are a consequence of strong correlations. We also show 
that the shift of the PL line is directly related with its 
asymmetry by a simple sum rule. 

Using Fermi's golden rule, the PL intensity at zero 
temperature is given by 
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where to is the PL frequency, |z) is the initial (trion) state, 
|/) denotes the possible final states, Ej is the energy 
of the state \j), and the relevant part of light-matter 
interaction can be written as 
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where A is proportional to the vector potential, d\ cre- 
ates an electron at the QD (with spin a), and p a an- 
nihilates a valence electron at the QD, with the same 
symmetry as that of the corresponding absorbed photon. 

The experimental results indicate that the hybridiza- 
tion of the trion state with the continuum is small. In 
any case the hybridization between holes and extended 
states favors a singlet state. Taking the limit of vanish- 
ing hybridization of the VWF of Varma and Yafet 15] for 
this case [in analogy to Eqs. (JSj) , (J9j) and the discussion 
below them] , one obtains [19( 
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structure in either Vk or p(oj) has to be assumed [13|. In- 
stead, as we show below, the three main experimental 
results mentioned above appear naturally when the in- 
teractions are properly taken into account. 



We solve the problem using many-body VWF's [15j, 
Il8j |. Similar VWF's were used successfully to study the 
interaction between spins of two QD's [21(. We extend 
previous approaches to the case of finite U. To simplify 
the discussion, we begin assuming U — > oo, which is a 
good approximation when E d is near the Fermi energy 
£f = £fc F , which we take as zero in what follows. For 
U — > oo, the VWF proposed by Varma and Yafet is a 
good approximation for the correlated ground state \g) 
of H El 



where c£ creates a conduction electron in a spherical 
wave with wave vector k, kp refers to the Fermi wave 
vector and \S) is the many-body state with all valence 
states of the dot and all conduction states below the 
Fermi energy occupied. Replacing Eq. © in Eq. H}, 
the intensity becomes 



\g)=a\F} + J2^4ck,\F) 



(8) 



where \F) = c£ *c£ JiS}. Assuming constant Vj. and 
p(oj), the resulting approximate ground state energy E g 
comes from the solution of the equation 



/M = T |A| 2 pi(^-M, (4) 

where 

p 1 (w)=Y / \{fn)?S(^-E f ) (5) 

/ 

is the spectral density of the many-body state 

\l) = (dW kFl -d\cl Fl )\S)/V2, (6) 

obtained after the application of Hlm [Eq.©] on \i) 
[Eq.©]. 

Since in the possible final states |/), the holes are ab- 
sent and do not play any role, our task reduces to calcu- 
late pi for the Hamiltonian of the Anderson model, which 
reads 

H = E d J2dU<r + Ud\d 1 dld l + Y,^kcl a c k(J 

+ 5^(V fc 4cfc . + H.c.), (7) 

ka 

where E d = E° d - eV g /X (with A ~ 7 [H) is controlled 
by the gate voltage V g . 

It is easy to see that in the non-interacing case (U = 0), 
pi is given essentially by the spectral density of d elec- 
trons, which is a symmetric Lorentzian for the usual as- 
sumptions of constant Vk and constant density of ex- 
tended states per spin p{ui) [2(|. Therefore, to explain 
the observed behavior within a one-electron picture, some 



E d - 5 + (2A/tt) In [(W + 5)/5] = 0, (9) 

where 6 = E d — E g > 0, A = irpV 2 is the resonant level 
width and — W is the bottom of the band of extended 
states [22|. Note that in the limit V — * when E d < tp 
Eq.® gives \g) = Therefore, the state |1) can be 
regarded as the limit of the Kondo ground state, when 
the Kondo temperature Tk ~ exp(7ri? ( i/2A) — ► 0. 

In the thermodynamic limit, \g) is orthogonal to the 
state |1). However, it is easy to see that excited states 
can be obtained making electron-hole excitations on \g), 
which have energy near (l|ff|l) = E d [22] and hybridize 
with |1). This hybridization leads to a broadening of 
Px{oS), in a way analogous to a resonant level model, in 
which a localized state is mixed with a continuum [20(. 
Specifically we consider the states 

|fcO = a fc V>+£/3£W>, (10) 

k 

where 

\Vkk>) = -j=^2dlc ktr \e k <). (11) 

Here the subscript k (k') denotes extended states below 
(above) the Fermi energy. It is easy to see that opti- 
mization of the coefficients for each k' (minimizing the 
energy), leads to a' — a, f3 k = /3 k , independently of 
k'. In addition, |(l|#|fc')| 2 = |al^| 2 = (I - n d )\V k ,\ 2 , 
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FIG. 1: density of the state |1) as a function of frequency for 
several values of Ed- For the lower values of Ed the amplitude 
was reduced by the factors indicated in the figure. 



where rid = (g\ d^d^g) is the electronic occupation 
at the QD. From the effective resonant level model, one 
obtains a broadening of the PL line given approximately 
by w = (1 — nd)A. This is a simple, elegant result which 
agrees qualitatively with experiment near the interme- 
diate valence regime (IVR) \E& — cp\ ~ A, and pre- 
dicts an exponentially small width in the Kondo regime 
€f — Ed, Ed + U ~ £f 3> A. However, since double oc- 
cupancy is neglected, this result is clearly wrong in the 
other IVR \Ed + U — ep\ ~ A, while in the Kondo regime, 
virtual fluctuations through double occupied states were 
neglected. 

To improve this calculation, wc first replace the states 
\ykk') m Eq. (fTPj) by others \ykk'}, which take into ac- 
count double occupancy 
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where the coefficients arc determined minimizing 
(ykk' \H \ykk')i f° r fixed k and k' . Second, in a similar way 
to the resulting states \k'), we construct linear combina- 
tions | As) appropriate for the other IVR \E d + U— tp\ ~ A. 
This is done easily using a symmetry property of H: 
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with 



T 

efc/ = — efc (for e_F = 0), which changes -E d into E d + U 
and therefore interchanges IVR's. For example using this 
transformation |efc/) — > djdj J2<y c k F a Ckcr \S) / V^- Third, 
since the states \k') and \k) are not orthogonal, we dis- 
cretize them and use a symmetric Lowdin procedure to 
orthogonalize them [23}. Finally, including the state |1), 
we numerically diagonalize H and calculate the spectral 
density using Eq. ([5|). 

The parameters were chosen as follows: A is our unit 
of energy. It should de of the order of 1 meV to give the 



FIG. 2: Half width at half maximum of the density of the 
state |1) as a function of Ed- 



observed PL width for Ed ~ £f- Guided from experiment 
[l3j |. we took U — 24A. Finally we took a broad flat band 
symmetrically placed around tp of width 2W — 100A. 
The results depend only logarithmically on W [see for 
instance Eq. 0] and the band edges are very far form 
the region of interest. Therefore there are no particular 
features of the band or hybridization that can affect the 
results. To discretize the states \k') and \k), the band was 
divided in at least 2000 points with decreasing energy 
separation towards the Fermi energy [2~j | . 

In Fig. [T]we represent pi(w) as a function of frequency 
lo for various values of Ed- This density with the fre- 
quency axis inverted and shifted by the energy of the ini- 
tial state Ei — Es [22} corresponds to the PL spectrum 
[see Eq. 0]. Due to our choice of a symmetric band, 
and the above mentioned symmetry property of H, the 
densities for Ed and — E d — U are identical. Therefore, 
one can restrict the study to Ed > —U/2 = — 12A. Note 
that in most curves an asymmetry with a high energy tail 
(corresponding to a low energy tail in the PL) is evident, 
in agreement with experiment. We return to this point 
below. For Ed — or slightly negative, there are two rel- 
ative maxima in pi(w). This structure is already present 
using the simplest VWF's given by Eqs. (|10p and pip . 
The peak at smallest energies corresponds to excitations 
near the ground state. At low energies, one expects that 
Px{oS) has a power-law behavior characteristic of x-ray 
edge singularities 14 , [2^| which is not captured by our 
variational approach. This structure might be washed by 
a temperature or other broadening effects not included 
in our calculations. This is likely the reason why a struc- 
ture with two peaks has not been so far reported. Note 
that while the density near the ground-state energy can 
be described more accurately with the numerical renor- 
malization group 1J, |25[ , this technique does not resolve 
accurately the structure at higher energies [25|, [27| ■ 

One striking feature observed in Fig. [T] is the dra- 
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FIG. 3: Shift of the maximum of the density of the state |1) 
as a function of Ed- 

matic narrowing of the peak as the system leaves the 
IVR and enters the Kondo regime, in Fig. [2] we rep- 
resent the half width at half maximum w as a function 
of the dot level, which is related to the gate voltage by 
E d = -eAVg/7. The slight increase of w for AV g = 
(and also E d = —24 A) might be related to the two-peak 
structure of pi(w) near these values of Ed and we believe 
it is not too significative. Instead, for the symmetric case 
Ed = — 12A, we obtain ro = 5x 10~ 4 A, corresponding to 
a decrease of the line width in three orders of magnitude. 
This is consistent with different experiments [H, [ll|, [l]| . 

In Fig. [3] we plot Em, the energy hoj for which 
Pi(oS) reaches its maximum value as a function of Ed = 
—eAVg/7. The shift Em — Ed is always negative (the 
peak is displaced towards the ground state energy). 
This means that the energy of the final state decreases 
and therefore, the position of the PL peak is always 
blueshifted with respect to the situation without hy- 
bridization. However, the non-monotonic behavior ob- 
served in Fig. [3] implies that as AV g increases starting 
from negative values, the PL peak first blueshifts, but 
after AV g = 0, the peak redshifts, in agreement with 
experiment [l3| . The fact that the maximum shift is ob- 
tained near the IVR's might be expected from the fact 
the maximum energy gain in the ground state is obtained 
at those regimes. 

It is easy to see from Eq. that J d(hui)pi(u>) = 
(l\H\l) = Ed [22]. In words, the center of gravity of 
the density states remains at E d when the QD electron 
hybridizes, while its maximum displaces at a lower energy 
Em- Therefore, to satisfy the sum rule, pi(w) should be 
asymmetric with a high energy tail, and the degree of 
asymmetry is directly proportional to the magnitude of 
the shift \E M - E d \. 

In summary, using an impurity Anderson model to de- 
scribe the final states of a photoluminiscence experiment 
in a quantum dot, with the simplest possible assumptions 
for the density of extended states and their hybridization 
with the dot electrons, we obtain a blueshift and cor- 



related asymmetry of the line shape, and a line width, 
which explains the main features of recent experiments 
in the whole range of applied gate voltage. 
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